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The problem of the figure of the earth has, as far as I know, only 
received one solution, namely, that of Laplace/^' His solution involves 
an hypothesis as to the law of compressibility of the matter forming 
the planet, and a solution involving another law of compressibility 
seems of some interest, even although the results are not perhaps so 
conformable to the observed facts with regard to the earth as those of 
Laplace, t 

The solution offered below was arrived at by an inverse method, 
namely, by the assumption of a form for the law of the internal 
density of the planet, and the subsequent determination of the law of 
compressibility. One case of the solution gives us constant compres- 
sibility, and another gives the case where the modulus of compressibility 
varies as the density, as with gas. 

It would be easy to fabricate any nuinber of distributions of 
density, any one of which would lead to a law of compressibility 
equally probable with that of Laplace ; but the solution of Clairaut's 
equation for the ellipticity of the internal strata of equal density 
seems in most cases very difficult. Indeed, it is probable that Laplace 
formulated his law because it made the equation in question integrable,. 
and because it was not improbable from a physical point of view. 

The following notation will be adopted : — 

For an internal stratum of equal density let — 

r be the radius vector of any point, 

a the mean radius of the stratum, 

e the ellipticity, 

w the density, 

colatitude from the axis of rotation, 

p the hydrostatic pressure at the point r, 0. 

For the surface let r, a, ^, It) denote the similar things. 
Let M be the mass of the planet, p its mean density, to the angular 
velocity of rotation, m the ratio uo^j^irp. 

* Since this paper was presented I haye seen a reference to a paper by tlie late 
M. Edouard Eoche, in vol. i of the Memoirs of the Academy of Montpellier (1848), 
in which the problem is solved, when the rate of increase of the density varies as 
the square of the radins. See Tisserand, " Comptes Eendns," 23rd April, 1883. 

t Laplace's hypothetical law of compressibility arises from a law of internal 
density for which the problem had previously been worked out, as an example, by 
Legendre. See Todhimter's "History of the Figure of the Earth," vol. ii, pj). 117 
and 337. 
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Let h be tlie ratio of tlie density of the stratum a to tlie mean 
density of all tlie matter situated inside tliat stratum, and 6 the 
surface value of h.^ 

Let 0, A be the greatest and least principal moments of inertia of 
the planet about axes through its centre of inertia. 

Let ce be the ellipticity which the surface would have if the planet 
were homogeneous with density p, so that CB-=:^m. 

The condition that the surface of the planet is a level surface is 
satisfied by — 

0-A=-|Ma3(c~im) ,..,.. (1). 

The condition that the internal surfaces are also surfaces of equi- 
librium demands that e should satisfy Olairaut's equation— 



( 



^^^ 6A^PWc^(* + 2^(;a3(^ + ?-lr=0 .... (2), 



dcb^ a^ / J \da a. 



It may be proved from (2), and the consideration that iv must 
diminish as a increases, that e cannot have a maximum or minimum 
value. 

Also it may be shown that the constants introduced in the integral 
of this equation must be such that — 

ce 5m 1 <:? / ox /o\ 

__= — = (ea^) ....... (S). 

t 4e 2tada^ ^ ^ ^ 

when a is put equal to a after differentiation. 
The mean density is given by — 



aV=3 



a 
wa^da ., . . (4), 





And fe=:~. 



Neglecting the ellipticity of the strata, we have the moment of 
inertia about any diameter of the planet given by™ 



=:-f7ri waMa ....... (5).. 



The ratio of (1) to (5) gives the processional constant. The pres- 
sure and density are connected by the equation — 

I {~-^-\-4^7rwcb\da-\- ~\ wa^da:=:0 .... (6). 
J a\w da / a jo 

* Ic, It, are the reciprocals of/, f, according to the notation adopted in Thomson^ 
and Tait's ''Nat. Phil." (edit, of 1883), § 824. 
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Now if t»y be a function sucL. that wdp=.dm, the dilBPerentiation of 
(6) leads to — 



dw 4i7r 
, — -i-__ 

cla a^ , 
and a second differentiation to — 



/* 



a 



iva^da=0 ....... (7). 





daP' 



{ma)-^4i7rwa=:0 ....... (8). 



It is well known that Laplace assumed that the modulus of com- 
pressibility of rock varies as the square of the density. Since this 
modulus is wdpjdiv, Laplace's hypothesis makes m proportional to w, 
and the equation (8) is at once soluble. 

After the determination of to as a function of a, the solution of all 
the other equations follows. 

In this paper I propose to find a new solution, and to compare the 
results with those of Laplace. 

In order to simplify the analysis let the unit of length be equal 
to the mean radius a of the planet, and the unit of time be such that 
the surface density hi of the plant is also unity. 

Kow let us assume that the law of internal density is — 

Then the mean density of all the matter lying inside of the stratum 
a is cr^^jil—^n). Hence, by definition we have — 

Thus we see that ^ is a constant for all strata, and therefore also 
for the surface. In Laplace's theory h is variable. With our assumed 
law of density and the special units, p the mean density is equal to 
the reciprocal of k. 

It is clear that n must be positive, otherwise heavier strata lie 
above lighter, and it must be less than 3 in order to avoid infinite 
mass at the centre of the planet. 

N'ow let us find the law connecting pressure and density, and the 
modulus of compressibility. 

Equation (7) becomes 



,.^__4_— 



3" T r~> 



da l—-^fb 



sand by definition of w and the assumption (9), 



dp 4^ a 



1-S?2 



da 1—^n 
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Integrating this, witli tlie condition tlaat tlie pressure vanishes at 
tlie surface, we have. 



2 






\~n 



3(l-w)(l-^w) 
whence the modulus of compressibility is 



[l-_^-2 ^ ]^ 



%0-^~^=-^7r W ^ n 



-9r 



\~n 



The case of %=1 is interesting; it gives a constant modulus of 
compressibility equal to Stt, and the law of pressure |9=27r log^. 

If n be less than unity the compressibility, or reciprocal of the 
modulus, increases with the density, which is of course physically 
improbable. If n be greater than unity and less than "3, the com- 
pressibility becomes less the greater the density. The assumed law 
probably does not give such good results as those of Laplace, because 
the decrease of compressibility with increasing density is not 
sufficiently rapid. 

The range of %=3 to ^=1 gives the results which possess most 
physical interest. 

In comparing results with those of Laplace there will be occasion 
to express the modulus as a length. ; that is to say, we are to find the 
length of a column of unit section whose weight (referred to the 
surface gravity of the planet) is equal to the force specified in the 
modulus. 

How if g be gravity 



3^ 



^ 1 fifJY^^ 

Hence the modulus is ^w^ » =^altt X - ( ~ I ** , the units a, ^ being 

reintroduced to give the expression the proper dimensions. Now 
galu is a pressure, and therefore the length of the modulus is 



2»-l 



- 1 - I *" . Thus the surface matter has a length modulus equal to aln. 

Now let us find the ellipticity of the internal strata. 
Substituting for w from (9) in (2), we have, 

da?' da a 

If the solution be assumed of the form 6=^GaP, jB must satisfy 

i3(/3-l) +2(3-^)/3~2?^=0, 

whence ^= — (f--t^)=l= ^/{(f)^— -^(S— w)}. 
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'Now n(S-—9i) is a maximum wlien it is equal to (f)^, and therefore 
tlie square root can never become imaginary. 

From the sign of the last term in the equation for ^, it is clear that 
one of the values of /3 is negative. Hence to avoid infinite ellipticity 
at the centre, the c corresponding to the negative root must be zero. 
Hence the solution of Clairaut's equation (2) is 

The surface value of -(ea^) is clearly '?^~|•+ \/ {(fy—n(3—n)} .. 

* CvCv 

Thus from (3) we have 

^=i(«_i)+iy{(|)3-„(3-»)}. 



Then substituting for w from (9) in (5) 

C = 



TT 



3(5-^) 

And since M=f7r/)=f7r --, we have for the precessional constant 

1 — pi 

C-— A ^—'yi/^ 1 \ 
= ft—-^m). 

C S-n ^ ^ 

J^I^ow let us collect these results, and express them in terms of Jc 
instead of 7i. The solution is 

And the mass inside of any radius a is 777-0^^'^. 

oh 

4-Qk 



/<2-3/0 



4— 6A; 



dp 4, 10^-3 A; 



W-:^=^f7r- 



dtv -^ 3/<l-/i;) 

And when ^=f, Jp = 27^1og^(;, iy-i_=27r, tu=:^a~^. 

dw 

The length of the modulus at the surface is 1/3(1 — li) of the 
planet's radius. 

--i-p+fv/{ap-/<i~/4^ 

z 

2 + 3/^ ^ 

O-A 2-4-3A^. 1 . 
3i^ ^ ^ ^ 
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Any value from unity to an infinitely small value may be assigned 
to 7^;, tliat is to say, we may have any arrangement of density from 
homogeneity to infinitely small surface density, but if h be greater 
than -| the compressibility increases with the density, which is 
physically improbable. 

The infinite density and infinite pressure, which occur in this 
solution actually at the centre, may be avoided by imagining the 
centre occupied by a rigid ;;spherical homogeneous nucleus, of very 
small radius ^a^ and of density \lhi)aF'^^~^^\ 

We have to compare this solution with Laplace's. 

For this case k is not constant, and its surface value is 6. 

Let 3-=.-^ where /c is a constant, being the arbitrary constant intro- 

IC 

duced in this solution ; and let 9 be the surface value of S-.^ 
The solution is — 

6 sin-S- 
«y sin 

And the mass inside of any radius am --~a^^. 

dw 



3(1-^ cot J) 

The length of the modulus at the surface is 1/(1 — ^ cot ^) or 3fe0"~^ 
of the planet's radius. 

C=[l-~6(l-~fe)^-2]fMa3. 
0-~A_ 1 ,,1 . 

"o~--Ti:6(i~fe)^-2^^ ^'^^• 

The following table gives the numerical values of the solution, 
together with columns for comparison with the results of Laplace's 
theory, for various values of the ratio of surface to mean density. 

^ See Thomson and Tait*s '' Nat. Phil.," 1883, § 824. 
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The case of 1^=0 gives the planet infinite mass at the centre, and 
the values are onlj inserted in order to complete the series. 



"k. 


Length mod. 

at surface 

m term.s of 

K as unity. 


where 


m 


Laplace 


C-A 


Laplace 

C-A 


where 




cc = 


t^ 


t 


0{t-hn). 


G{t-~lm). 




I'O 


00 


0-000 


1-000 


I'OO 


1-667 


1-67 


■— oo 


•9 


3-333 


0-132 


1-066 


l-OJf 


1 


741 


1-71 


-4-667 


•8 


1-667 


0-293 


1-147 


1-09 


1- 


833 


1-77 


-1-333 


•7 


1-111 


0-488 


1-244 


1-15 


1 


952 


1-83 


! -0-067 


•6 


1-^1-2 


0-722 


1-361 


I'U 


2 


111 


1-90 


1 0-333 


•5 


1-^1-5 


1-000 


1-500 


1-27 


2" 


333 


1-98 


0-667 


•4 


1--1-8 


1-322 


1 -661 


I'SIp 


2 


667 


2-07 


i -889 


•3 


l-f-2-1 


1-688 


1 -844 


1-41 


3- 


222 


2-17 


i 1 -048 


•2 


1--2-4 


2-093 


2-047 


1'48 


4 


•333 


2-28 


1-167 


•1 


l-~-2-7 


2-532 


2-266 


1-56 


7-667 


2-40 


1-259 


•0 


l-v-3-0 


3-000 


2-500 


1-65 


00 


2-65 


! 1 -333 


•667 


1-000 


0-562 


1-281 


1-17 


2 -000 


1-85 


1 jpoc log w 


•333 


1-^2 -0 


1-562 


1 -781 


1-39 


3-000 


2-13 


1 1 -000 



Note. — The values in the two columns applicable to Laplace's theory were found 
by graphical interpolation from a series of values given in " Month. I^ot.," R.A.S.j 
Dec, 1876, or Thomson and Tait, "Nat. Phil." (1883), § 824'. 

In Laplace's theory poc {w^ — l), and the modulus of compressibility oc w^. Li 
the present theory the modulus oc w^. 

The value ^ = '667 corresponds to constant compressibility, and l: = -333 to 
gaseous compressibility. 

One of the grounds on which Laplace's solution is held to be satis- 
factory is that if we take the value of ce, as determined by the known 
angular velocity and mean density of the earth, and the value of z as 
determined by geodesy, and find the value of % the ratio of surface to 
mean density, which corresponds with the ratio cejz^ this same value 
of fe is found to give a proper value to the coefficient of z — ^m, so as 
to obtain the observed processional constant. To be more precise, m 
is found to be 1/289-66, which gives ci5 — 1/231'7, and z has been found 
to be approximately 1/295. These give ce/^=: 1-273, and this corre- 
sponds with fe=l/2-06 = *49. This value of % with the same values 
of z and w, gives the processional constant as '0033, and Leverrier 
and Serret give its value as '00327. 

E'ow it appears remarkable that almost as good a correspondence 
is obtainable from my solution. The value celz=l'27S corresponds 
with fe=:*675, and when fe = -675 the coefficient of z — ^m in the pro- 
cessional constant is 1-99, which gives the same processional constant 
•0033. 

This value of ^ corresponds very nearly with constant modulus of 
compressibility, and with pressure determined by jp=:27r log iv. 
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It is claimed in favour of the Laplacian liypotliesis tliat it corre- 
sponds to a surface density wMcTi is nearly a lialf of the mean 
density of the earth, and that we know that average rock has a den- 
sity of about 2*8. Also it is pointed out in Thomson and Tait's 
" Natural Philosophy " that the length modulus of compressibility 
of the surface rock is about 1/4*4 of the earth's radius, which is 
very nearly the observed length modulus of iron. 

These conditions are not well satisfied by the present solution, for 
the surface density is found to be 'B^S, or 1/1*48 of the mean density 
of the planet, whence the specific gravity at the surface is 3*7; and the 
length modulus at the surface is equal to the planet's radius. It is fco 
be admitted that this density is large, and that the substance is also 
highly incompressible. 

Thus in these respects the Laplacian hypothesis has the advantage. 
It seems to me, however, that too much stress should not be laid on 
these arguments. We know nothing of the materials of the earth,, 
excepting for a mile or two in thickness from the surface, hence it is 
not safe to argue confi-dently as to the degree of compressibility of 
the interior. There seems reason to believe that there is a deficiency 
in density under great mountain ranges, and this would agree with 
the hypothesis that our continents are a mere intumescence of the 
surface layers. 

According to this view we might expect to find a rather sudden 
change in density within a few miles of the surface. N'ow in any 
theory of the earth's density such a sudden change in the thin shell 
on the surface could not be taken into account, and the numerical 
value for the surface density should be taken from below the intu- 
mescent layer if it exists. Hence it is not unreasonable to say that a 
solution of the problem, which gives a higher surface density than 
that of rock, lies near the truth. I do not maintain that my solution 
is as likely as that of Laplace, but it is not to be condemned at once 
because it does not satisfy these conditions as to the density and com- 
pressibility of rock. 

The two cases which are given at the foot of the above table each 
possess an interest, the first of constant compressibility, because it 
corresponds with the case of the earth, and the second of modulus of 
compressibility varying as the density, because this is the gaseous 

law. 

With constant compressibility the internal ellipticity varies as the 
662 power of », or nearly as the square root of the radius; with 
gaseous compressibility it varies as the 1*562 power of a^ or nearly as 
the square root of the cube as the radius. 

A numerical comparison of the case of constant compressibility 
with Laplace's solution for fe=J gives the following results :— - 
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a = 








\n. 


i a. 


|a. 


a. 


(Laplace) — = 

(Constant e 
compress.) ^ ~ 


•812 
•459 


•844 
•693 


•902 
•851 


1-000 
1-000 



Thus tlie Laplacian solution attributes mucli higlier ellipticitj to 
tlie internal strata. The solution with constant compressibility iu 
fact gives so large a proportion of the mass in the central region, that 
attraction has a greater influence compared with rotation, than in the 
solution of Laplace. 

[P.S. — If, as is not improbable, the increase of density in the interior 
of the earth is due rather to the heavier materials falling down to the 
centre than to great pressure compressing the material until it has a 
high density, then the determination of a modulus of compressibility 
would be fallacious, and it would be more logical to leave the 
expressions for the pi'essure and the density both as functions of the 
radius, without proceeding to eliminate the radius and to form an 
expression for the modulus of compressibility. I owe this suggestion 
to a conversation with Sir William Thomson. — December 19^ 1883.] 



